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Abstract 

Functional conjugation methods are used to analyze the global structure of various renormalization group trajectories, and 
to gain insight into the interplay between continuous and discrete rescaling. With minimal assumptions, the methods produce 
continuous flows from step-scaling a functions, and lead to exact functional relations for the local flow j3 functions, whose 
solutions may have novel, exotic features, including multiple branches. As a result, fixed points of a are sometimes not true 
fixed points under continuous changes in scale, and zeroes of /? do not necessarily signal fixed points of the flow, but instead 
may only indicate turning points of the trajectories. 



I. INTRODUCTION 



The renormalization group (RG) of Gell-Mann and Low [lij], and of Stueckelberg and Petermann [2~il]. has an 
elegant mathematical expression in terms of the functional conjugation (FC) methods of Ernst Schroder [23|. This 
expression provides a powerful tool to describe the behavior of physical systems under either infinitesimal or finite, 
perhaps large, changes in scale. While this fact is often overlooked, and not usually invoked in the solution of various 
problems posed in the RG framework, it is readily apparent upon reading [lij (see especially Appendix B; also see 
|19j [) and surveying the literature on functional equations [161 ]. Moreover, it may be profitable to bear in mind the 
logical connections between these two subjects when considering the step-scaling approach in lattice gauge theory 
Pri20j. where the power and utility of the methods are manifest. 

In previous work we have discussed how dynamical systems, defined on a discrete lattice of time points, may 
be smoothly interpolated in time through the use of solutions to Schroder's celebrated functional equation. Here 
we discuss the same methods in the context of the renormalization group. We examine in detail the connections 
between differential (local) rescaling and finite (global) changes in scale. We interpolate various step-scaling functions 
to obtain trajectories under continuous change of scale, with emphasis on the consistency imposed by the analytic 
properties of couplings in the presence of UV and IR fixed points. From this point of view it is possible to obtain 
novel features for RG behavior. In particular, multi- valued Callan-Symanzik /3 functions [1, 0, El] are commonly 
encountered in the local RG flow equations, even when interpolating elementary, polynomial step-scaling functions, 
with interesting consequences involving fixed points, cycles, and even chaotic evolution under changes in scale. 

In section II we describe functional conjugation methods relevant to RG analysis and apply them to the study of 
selected trajectories. In section III we consider a physical illustration of fixed point behavior drawn from numerical 
studies of lattice gauge theory [lj . In section IV we illustrate elementary limit cycle behavior in a model obtained by 
an extension of the standard BCS Hamiltonian [13] ■ In section V we briefly explain how further novel, exotic features 
can arise from basic step-scaling behavior, in general. Finally, in section VI we exhibit such features, including 
multi-valued /3 functions, limit cycles, and chaotic trajectories, using toy models based on the logistic map. Two 
appendices provide some connections to our earlier work on dynamical systems, and a few algebraic details for the 
lattice example. 



II. METHODOLOGY 



A. The renormalization group: Step by step 



Let us suppose the change in the coupling u is given for a discrete change in length scale by 

u H> a (u) , (1) 

where a (u) is the "step-scaling" function 0, H3- Typically, u — g 2 /47r, where g is the gauge coupling, although it 
may be convenient to incorporate other numerical factors into u, or even to take other functions of g 2 , depending on 
the problem at hand. The standard interpretation is to regard a (u) as a discrete At sampling of a renormalization 
trajectory, u(t), whose continuous evolution under changes in the log of the length scale, t, has proceeded from an 
initial u = u (t)\ t=Q . This notation for the initial u (rather than uq, say) is not only more convenient to express the 
step-scaling function, e.g. as the mapping {jTJ) , but also to write many relations that hold both for the initial u as well 
as more generally for all u (t) . We will point out several such relations in the following. 

The trajectory is assumed to describe an abelian i-flow with group composition given by simple addition of t 
arguments. So, for example, a (u) = u(t)\ t=1 , a (<r (u)) = u(t)L_ 2 , c -1 (it) = u(t)\ t __-,, etc. The local flow 
equation in terms of t has both familiar and more recondite forms, 

on 

-= J 8(u) = (lnA)*(«)/*'(u) ) (2) 

where (3 is the so-called Callan-Symanzik function [!,[25[ (which appeared earlier in (24|, Eq(4.25), under the alias h% p ) 
and \t is the so-called Schroder function [23| (both of which appeared in [l4j, under the aliases if) and G, respectively), 
and where 1/lnA (usually taken to be ±1) sets the scale of t. It is also implicitly understood that the system is 
underlain by a i-translation covariance so that ^ holds not just for u = u (£)| t=0 but also for u (t), provided of course 
that the RHS is also modified by p (u) ->• ft (u (i)) and * (u) -> * (u (i)). 

While the Schroder function is less well-known in renormalization theory, it is immediately expressed in terms of j3 
from the definition in rewritten as 

dln^(u)/du = (ln\)//3(u) . (3) 
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Thus a definite integral gives the total change in ^ brought about by a finite change in the coupling, 

$ (u 2 ) = A^''i 2 * (m) . (4) 

On the other hand, the exponent here is just t 2 — t\, the total change in t as u% — > u 2 , as follows from the first equality 
in ([2]). Therefore another way to express (j4|) is in terms of the evolution of the Schroder function under the flow of 
the coupling, u — > u (t) , 

^(u(t)) = A* *(it) . (5) 

This last relation reveals the fundamental role played by "J, and its inverse function $"\in the construction of 
trajectories for arbitrary changes in t. It follows from ([5]) that such global flow is given by [14[ 

u{t) = ^~ 1 (A t *(?i)) , (6) 

where 1/lnA sets the scale of t. The RHS of ^ is immediately recognized as just a change of variable, effected 
through a functional conjugation 23]. For us, in fact, the expression "Schroder functional method" is just a metonymy 
for functional conjugation. 

Indeed, the expression ^ is perhaps the most succinct way to appreciate that renormalization relates self-similar 
structures at different scales, inasmuch as the RHS is just a functional similarity transformation: o A* o f , 

Moreover, ([5]) shows that fixed points or limit cycles can arise in a model for real A if and only if "J -1 either becomes 
constant or else exhibits periodic behavior, respectively. 

The structure of © also makes the abelian t-flow of the renormalization group manifest, and it gives a formula for 
the step-scaling function, or any of its functional compositions, in terms of ^S. For example, for a (u) = u (t)\ t=1 , we 
have from ([5]) 

\y(u) = y(a(u)) , (7) 

a form known as "Schroder's functional equation with eigenvalue A." Presented in this form, for a given a(u), the 
problem is often to determine all allowed A and to find all solutions of the functional equation [2^ | . 
Alternatively, we may write ([7]) as 

cr(u) = *- 1 (A*(m)) . (8) 

In this form, the equation determines the step-scaling function in terms of 'J. In fact, it is useful to think of u (t) in 
([6|) as at (it) j that is, as a t-th continuous functional composition of er. For example, a 2 (u) = a (a (u)) as before, but 
now generalized to a (u) = cri/ 2 (o~i/ 2 (u)), etc. More generally, c^+ta (u) = a tl (cr t2 (u)) — just the expected RG 
abelian composition rule. 

At this point it is natural to ask, what is a simple physical model whereby ^ (t) — A*^? Well, dln^ (t) jdt = In A, 
so clearly In ^ is the variable of choice. Then the question becomes, for what model is the change in the coupling 
with scale a constant? An obvious answer is, the one-loop approximation for evolution of an inverted coupling, 1/g 2 . 
That is to say, if 

^0(i)=/Wop = ic 3 (t) , (9) 



then 

d 



dt Vfl 2 (*) 



= -c ■ (10) 



So, the physical interpretation of the Schroder function is clear: The log of \P is just the change of variable needed 
to convert the renormalization group flow for u into a one- loop flow for a re-defined coupling constant l/fl 2 . Thus, 

d ln*( U (t))=lnA ±(J—)=-c. (11) 



rfi dt \q 2 (t) 

The role of In ^ is put into deeper perspective through the following formal sequence of steps that evoke the method 
of characteristics for the RG. Making explicit the dependence of the trajectory on the initial u = u (£)lt=o as weu as 
on t, and making use of 

4- = P(u) 4- = On A) ni , ( 12 ) 
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along the trajectory, we have 



u(t,u)= u(t,u) = e t0M ^ u = e t(lnA) ^w u . (13) 



T=0 



But now, u = * 1 (* (u)) = * 1 (exp (In "J (u))), so the last expression reduces to a mere translation of the variable 
In (*(«)), 

e t(lnA)o (exp (In* (u))) = (exp (t In A + In * («))) = * _1 (A** («)) . (14) 

Thus (j6]) is recovered. These formal steps can be made precise by examination of T, the radius of convergence of the 
t series, in particular by a determination of the dependence of T on the initial u. How this goes will be illustrated in 
the examples to follow. 

There is some additional, subtle mathematical structure to take into account here, especially if we have in hand a 
series expansion for a: 

a{u) = au + (u 2 ) . (15) 

For example, if 4rU (t) = cu 2 (t), as is true for lowest order perturbation theory with u oc g 2 , then the exact solution 
for the trajectory is 

^ = I^< (16) 



where again on the RHS u = u (t)\ t=Q . In this case we have a = 1 in (I15[) . 

But this leads to A = 1, the well-known singular situation for the eigenvalue in Schroder's functional equation. By 
singular we mean that an attempt to solve ([7]) by Taylor series expansion about u — will fail, in general, when 
A = 1. Indeed, it is immediately verified that, for a = 1 in (|15[) . a nontrivial solution of ([7]) can not be found with 
* (u) given by a series about u = 0. 

This is easily circumvented, however. Instead of Taylor series about u — 0, all the relevant solutions have an essential 
singularity at u — 0, and in fact have Taylor series about u = 00. Explicitly, with A = exp (nc) in @, we see that 
*« (it) = exp (— is a family of Schroder functions for one-loop evolution, with arbitrary k: *„ (u (t)) = e KCt ^ K (u). 
So, for t — 1 and a (u) = u (t)\ t=1 , the eigenvalue for each * K solution is indeed A, and not the naive value 1 (excluding 
the trivial and uninteresting case where \l/ K =o = 1)- Defining a discrete step for another value of t simply rescales k. 

We will say more about the general structure of the functional approach to the RG, and the novel features that it 
has the power to reveal, in sections IV and V of the paper. But first, we consider: 



B. The p functional equation, with one- and two-loop examples 

What is the functional equation obeyed by the local /3 function? It follows simply enough from ([2]), or else from 
the definition in ([2} combined with (JJJ- Thus, 



(17) 



That is to say (cf. Eq (6) in ^ 



Ha{u)) = ^M Hu) . (is) 
du 



This is just the flow of a (u) along the characteristics defined by (fl~2j) . Note that all explicit reference to the eigenvalue 
A has dropped out of this equation, although it is still possible for A dependence to be induced through implicit A 
dependence in a (u), and therefore A dependence is implicitly understood for (3 as well. Also note that (TT8|) alone does 
not determine the overall normalization of /3. This normalization is determined by @ and ([7]). and it also introduces 
A dependence. 

As an example, again take (3 (u) = cu 2 , the one-loop result. For t = 1, we have a (u) = ± ™ , hence 

^P- = -^— 2 • (19) 
du (l-cu) 2 
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But we also have 



[3 (a (u)) = co~ («) = , = , /? («) . (20) 

(1 — cu) (1 — Clt) 

So, (|18p holds — the most important point to take away from this example being that a Taylor series solution about 
u = can work for the functional equation obeyed by ft (u), even though it does not work for '5 (it). It is of course 

redundant to do so, but we check this, for j3 (u) = au + ju 2 + bu 3 . Then a (u) — l ^ cu gives /3 (<r (w)) — ( da }^ ^ (u) = 

^ e ^" 1 ^j (— 6u 2 — acw + a). For this to vanish, it is necessary and sufficient that both a = and 6 = 0, while 7 is 
undetermined. 



1. Two-loop infrared fixed point 

For another, more interesting example, consider the two-loop perturbative approximation to /3 for a model with 
trivial UV and nontrivial IR fixed points. This example nicely illustrates how the normalization of (3 is determined 
in the functional approach. We may sweep various model dependent factors into the definition of the coupling, g, 
and the scale of t to write: 



Changing variable to 



211) becomes 



with solution \li 



y = 4 - 1 - C 22 ) 
g 2 



% = TT~> (23) 
at l + y 



y (t) = LambertW (yoe^ "') , (24) 



where LambertW is the inverse function for xexpa;. This is manifestly in FC form ([6]), with A = I/e and "J/ (y) = 
j/expy, where \& satisfies the functional equation 

1 V (y) = * (LambertW [- yexp(y)]) . (25) 



e \ \e 

From this, we immediately read-off the step-scaling function in terms of y. Switching back to the original coupling 
g the solution (|24l) gives 

9 2 (*) = 77 1 x , N • (26) 

I + LambertW 

A typical trajectory is shown here. 



((i-l)e-^) 
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-10 



10 

t 



A 2-loop trajectory with g (0) = 1/2. 

In this example, ^ can be obtained by Taylor series solution of the functional equation, Q25|) . only it is a Taylor 
series about the non-trivial fixed point at g = 1, i.e. y = 0. Nonetheless, for this two-loop example, one may forego 
the series solution of (|25j) and just solve it by inspection, upon noting that LambertW -1 (z) — ze z . 

For this same example, the /3 function can be obtained from the step-scaling function by series solution about either 
fixed point, g = or g = 1. In terms of the variable y with tr (y) = y (£)| t=1 , we have 

LambertW (ye^ 1 ^) 



a (y) = LambertW y — yexp (y) 
This leads to a typical plot of a near the fixed point. 



do- (y) 

dy 



1 



y (1 + LambertW {ye~ 1+ y)) 



(27) 




Two- loop IR fixed point exhibited by a (y) (solid red) and a 1 (y) (dashed red) versus y. Light gray curves are the 

functional square roots of a and er -1 . 

Of course, with all the additional information about the actual trajectory implicitly built into this closed-form 
expression for a, we may also forego an actual series solution of (|18[) and solve it too by inspection. Explicitly writing 
out the functional equation as 

LambertW (ye~ 1+y ) l + y 



P (LambertW {ye~ 1+y )) 



(28) 



(1 + LambertW (ye~ 1+ y)) y 
a solution is obvious, namely, /3 (y) oc jq^y- The constant of proportionality is then given by In A = In (1/e) = —1, as 
in ((2|), and thus ([23)) is recovered. 



G 



2. Mobius transformation form 

This is a simple Pade approximant (l2j , sometimes known as an "exact" j3 function [22j , and is but a minor variation 
on the previous two-loop example. For constants a, /?, 7, and 5, consider 



dg = 1 3 a + /V 
alt 2 ff 7 + <5 5 2 



5 d — • ( 29 ) 

at z 

Upon changing variables to 



ad — f3j ' 

the equation becomes just like (|2"5|) with only a change of scale for t: 

dy _ -a 2 y 
dt aS — /?7 1 + y 

Thus the solution is 



7(/W, 2 ) 



(31) 



?/ (i) = LambertW (yoe yo Q<5 -^ 



This is again of FC form with eigenvalue A = e aS -^ . In terms of the original variable the RG trajectory is given by 

9 2 (*) = 7 7 — 7 7 \ \T • (32) 

1 - (S^ftr) LambertW ^ + « j exp ^ + ° J - ^^)) 

5. One-loop geometromorphosis 

This describes the renormalization flow of geometry from a fiat manifold in the UV towards a fixed, nontrivial 
manifold in the IR. The trajectories describe evolution of matrices, 

d9ab (t) - P ab (*) , 9ab (t) = 9ab (0) + f dr p ab (r) . (33) 



dt j 

For RG flow of this type, we would expect a tensor version of the Schroder equation to be operative, 

* o b(g(l))=A ac * cfe (g(0)) , (34) 

where g a b could include the torsion potential as well as the metric (in which case g a b 7^ <7ba) and where the step-scaling 
function is now a RG transported g a b, 



g ab {l) = A ac g cb (0) , A QC = 8 ac + (J^ dr f3 ad (t) \ g db x (0) 



(35) 



To simplify the discussion, and to be explicit, consider the 3-sphere er-model with torsion Q, with © proportional 
to the square of the radius of S3. In this case, the one-loop renormalization of the metric boils down to just a change 
of & with length scale: 

d& _ 1 / 1 
~dt ~ 2 \& 



= o(=5- 1 ) ' ( 36 ) 



The solution of this one-loop evolution equation is given implicitly by 

© (*) + 1 „-26(t) _ Jt & + 1 



e (t) - 1 6-1 

That is to say, ([5]) and ([7]) have eigenvalue A = e, with explicit Schroder function 

6 + 1 



* (e) = t^T e " 2S ' (38) 
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while \fr 1 is only implicit. Note that ^ (6) > when (3 > 1. The 3-sphere squared-radius RG evolution is then 
given in FC form by 

&(t) = -if- 1 (e'*(6)) , (39) 

This defines implicitly the step-scaling function a (©), as 6 (1), say, shown here. The fixed point a (©*) = 6* = 1 
is centered in the small black circle. 




One-loop "geometrostatic" IR fixed point exhibited by a (6) (solid red) and a 1 (©) (dashed red) versus ©, and a 

few other fixed t slices of the at (©) surface (thin gray curves). 

The implicit function has no name, as far as we can tell, although it might be classified as a generalization 
of the Lambert function. In any case, we may construct the inverse function "f^ 1 (z) through the usual graphical 
methods, to find two branches. 

3 T 




Two branches of ^ 1 (z), in green, and various approximations (orange dashes). 
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The UV fixed point — a 3-sphere of infinite radius — corresponds to the verticai asymptote of the upper branch of 
$ , while the IR fixed point — a 3-sphere of unit radius — corresponds to the horizontal asymptote (thin light 
gray line in the Figure) of that upper branch. Also as shown in the Figure, it is straightforward to construct series 
approximations around the point (z, (z)j = (0, —1), for the lower branch, 

* Wr branch (*) = "1 " 2e" 2 * + 6e"V - 26e~ e z 3 + ^197 e - 8 z 4 - 722e- 10 z 5 + O (z 6 ) , (40) 
as well as around the point(s) (z, ^ r_1 (z)j — (±00, +1), for either the upper or the lower branches, 

« = 1 +2 «" (1) - (;)' + ^" (;)' - 1 1 ^-" (;)' + (;)' + ((;)*) ■ <«» 

The second of these series follows from the first through the transformation (\f r , ©) — > (A, — ©) applied to the equation 
for These series results for are representative of situations where simple, closed- form expressions are not 
readily available. A similar situation often arises when we have: 



III. A LATTICE GAUGE THEORY MODEL 



In this section we consider a physical illustration of fixed point behavior drawn from numerical studies of lattice 
gauge theory The j3 and a functions in question are those for a non-abelian gauge theory with 12 flavors of 

su (3) color triplets. 



A. Parameterizations 



It is convenient in lattice gauge theory to approximate the physical coupling at length scale £, u = 1/g 2 (£), 
parametrically in terms of the bare lattice coupling, s = l/<?o> as a series [l|, 



The step-scaling procedure [1, [2(| then gives the coupling at length I as a function of u. Of course, this is also a 
function of the bare lattice coupling parameter, s, but with a different series expansion, in general, perhaps even with 
a different order for the series (depending on the choices made in the numerical computations). 



«(«> = — = ' ! ' - >.'.< [ >-l 

Note that these series are arranged to have a common zero-coupling limit, as s — > 00, with both (l4"2l and (|4U)) becoming 
the identity map in that limit. For this parameterization, the functional equation (|18l) for j3 (u) may be written as 

a 1 1 1 \w du(s) da(u(s)) . . 

P(^(u(s))) = /3(u(s)) -j^ 1 , ( 44 ) 

where 

du(s) „ v"-* , . -\ //a 1 1 da(u(s)) \— \ , . , r , l 

-^ i=i +Eo- i )^w7 and y ^ =1+^0-1)^-^)-. (45) 

It is sensible from the stand-point of perturbation theory to consider a similar expansion for j3. Thus we write 
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Using this series along with (|42j) . (|43|) . and (1451) . and expanding both sides of (|44j) in powers of 1/s, we obtain recursion 
relations for the coefficient ratios b n /bo. The O (1) and 0(l/s) terms on LHS and RHS of (|44| match identically, 
but the O terms for n > 2 give expressions for b n -i/bo in terms of the Cfc<„. For example, we find 

h = c 2 (L) - c 2 (£) 62 _ 2c 3 (L) - 2c 3 {£) - (ci (L) + ci (I)) (c 2 (i) - c 2 (f)) ^ 

6 Cl (L) - a (£) ' 6 

etc. The overall normalization of /3 is not determined by (|44|l . of course, and in fact there is no information in the 
expansions (j4"2"|) , (|4"3"|) , and ([4"5"j) that allows determination of bo by Taylor expanding about s = 00. (This is related to 
the essential singularity in the Schroder function solutions at zero coupling, as mentioned earlier, (fTSj) et seq.) Rather, 
we must fix 6 by other considerations, the obvious choice being to use perturbation theory. Another possibility is to 
use lattice data and expand, not about zero coupling, but about a nontrivial fixed point, if available. The expansions 
are similar to those we have just given, and are collected together in the Appendix, §IX. 



ci (L) - ci U) 



B. Numerics 



Considerable effort is needed to properly take the continuum limit where the lattice spacing a goes to zero. However, 
for purposes of illustration of the various functional methods described here, we will not concern ourselves with those 
complications. Rather, we will simply take some of the raw numerical data in [l| for the expansion coefficients 
appearing in (|42p and (|43|) . and note with amusement that such a naive ab initio computation of the j3 function 
matches very well with two- loop perturbation theory, upon overlapping the two results. 

We choose to consider the L = 8a and L — 16a data from [l| for 12 flavors of su (3) color triplets, a model widely 
believed to have a nontrivial IR fixed point near g 2 — 5. We take the data at face value, without regard for any 
statistical or systematic errors, and we use this data for the parametric definitions of u and a (u). 

u (s) = s [ 1 - 0.4092 [ - ) + 0.192 



SI \s 




2 

s I \s I \s 



a (u (s)) = s I 1 - 0.467 I - ) + 0.154 [- \ - 0.164 I - I I . (48) 



For illustration purposes, this will suffice; but, hopefully, the procedures we follow will be useful in future, realistic 
lattice studies. 

In any case, we plot o~(u), parametrically, versus u to display the c(u*) = = 1/g 2 fixed point encoded in the 
data, as given numerically by: 

9* = tt^t = 5. 56 at s = 0.739 , i.e. gl = = 1. 35 . (49) 

y* 018(J yo* 0J39 v 1 

The two series (|48j) were constructed to become parallel curves in the zero coupling limit, as s or u — > 00. However, 
the curves have significantly different approaches to the nontrivial fixed point, hence its existence, as is evident in the 
Figure. The fixed point u* = 0.180 is centered in the small black circle. 
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a (u) (solid red) and a 1 (u) (dashed red) versus u, and a few other fixed t slices of the at (u) surface (thin gray 

curves) . 



Now, if we use the s-parameterization formalism discussed above, we obtain for example: 

-i = 0.657 , = -20.2 . 
bo b 



(50) 



This certainly looks useless (i.e. like a divergent series for g 2 <J 0.03) and is very much in disagreement with 
perturbation theory, as given to 2-loops for the 12 flavor su (3) model by 



d 

—u = -2 
dt 



50 1 



(4tt) z (4tt)' 



(51) 



such that bi/bo 



-0.106. But if nothing else, (|50p confirms that the data in (|48l) certainly did not come 



just from transcribing perturbation theory. 

On the other hand, if we reparameterize the data around the IR fixed point, defining 

- = (1 + r) gl , 

and expand various quantitites to and including O (r 4 ) , as described in detail in an Appendix, we find 

u (r) = 0.180 - 1. 52r - 1. 258r 2 - 3. 25r 3 - 0.408r 4 - 0.739r 5 + O (r 6 ) , 
a {u (r)) = 0.180 - 1. 13r + 0.439r 2 - 0.739r 3 + 0.739r 4 + O (r 5 ) . 



(52) 



(53) 



Admittedly, the second of these r-series is less firm, since we only have a [u (s)) to O ((i) 3 ), although the diffe rence 

will turn out to be slight over the region where the series is reliable. We describe below the changes encountered 
from truncating the a (u (r)) series at O (r 3 ) versus O (r 4 )- 

Using again the functional formalism to determine the j3 function from a, only now as implemented in terms of the 
expansion about the nontrivial fixed point (for details, see the Appendix), we find the more sensible expression 



P{u) = p (u* + w) = bi x (w - 3.14w 2 + 1. 86w 3 + 4. 02w 4 ) , where bi = In (0.745) = -0.294 . 



(54) 
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This is the best one can do given only the data in (|48|) . So, small, positive initial w will decrease as t increases. 
That is to say, u> u* will decrease down to it* as t increases. Or, since u — 1/g 2 , g 2 < g 2 will increase up to g 2 as t 
increases. Increasing t corresponds to increasing length scale, so the fixed point is indeed an IR one. This is clear 
from computing the trajectories that follow from (j54")l . Here is an example of g 2 (t), with g 2 (0) = 3. 




-20 -10 10 20 30 

t 



Lattice (solid blue) and 2-loop (dashed blue) g 2 (t) trajectories, with g 2 (0) = 3. 

In terms of both slope and curvature (cf. cubic splines), the trajectory joins very smoothly with that obtained from 
two-loop perturbation theory for 12 flavors of quark color triplets. 

Further comparison can be made to trajectories determined by other approximations. It is more direct, however, 
to just compare /3 functions, now that we have in hand (|54[) . In addition to the two-loop approximation, which is 
unchanged by analytic redefinitions of the coupling, we also compare to three-loop (minimally subtracted, as well as 
that based on the Schrodinger — not Schroder — functional) and to four-loop results from perturbation theory (see 
[l| for the literature dealing with these approximations). Note how well the two-loop result joins smoothly to our 
naive, but ab initio determination of f3 using the lattice data and functional methods, at g 2 — 3. We stress that there 
have been no adjustments in the normalizations to facilitate this match-up. The other, higher-loop approximations 
do not fare nearly as well in terms of matching-up with the lattice /3, although they do give estimates of g 2 more or 
less in line with (14T)1) . 



beta 




Various /3 (g 2 ) functions: 2-loop (black), 3-loop MS (green), 3-loop SFR (orange), 4-loop MS (blue), lattice (red). 

Finally, we contrast the differences resulting from truncating the o~(u(r)) series at O (r 3 ) versus 0(r 4 ). The 
resulting changes in (g 2 ) , as obtained from the step-scaling function, are shown below. The difference is slight near 
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g 2 = 3, but becomes a sizeable disagreement for g 2 iS 2.5, at which point the disparity is comparable to that between 
(|54|) and the higher-loop approximations at g 2 = 3. In fact, the O (r 3 ) truncation even gives another, spurious(!) 
nontrivial fixed point at g 2 sa 1.7. Simply put, for g 2 is 2.5 the O (w 4 ) series in (|54| is woefully inadequate. A 
better approximation to the solution of the /3 functional equation must be used for smaller g 2 . 



beta 




Two-loop (black) and lattice (red) /3s, the latter to both cubic (dashed) and quartic order. 

While some of these numerical coincidences may very well be little more than artefacts of the data selected, it is of 
interest to see if this match-up between perturbative results, and those obtained for j3 by functional methods, persists 
given more accurate lattice data and a thorough analysis of numerical errors. The above discussion is an illustrative 
application of the functional methods introduced here, rather than an endorsement of specific results. 



IV. A MODEL WITH A LIMIT CYCLE 



In this section we illustrate elementary limit cycle behavior in an exactly solvable lattice model, albeit one- 
dimensional, as obtained in pjj by an extension of the standard BCS Hamiltonian to include a term that breaks 
time-reversal invariance. A model with similar RG structure was studied in 15]. As before, we emphasize the 
connections between discrete and continuous scaling, and the underlying functional relationships. 

The physics of the model is explained in the extensive work of LeClair et al. and only the key mathematical features 
will only be summarized here. The Hamiltonian for the model is 

{(g + ih) e if ej > e k 
ge if ej = e k , (55) 
(g - ih) e if Sj < e k 

where bj and bj denote the usual Cooper-pair annihilation and creation operators, and where e = \ (sj+i — £j) is the 
single-particle level spacing. 

For a large number of system sites, the renormalization of the dimensionless couplings g and h under a change in 
system size L is given by 

d 9 2 , ,2 



g + h , h = constant , (56) 



d\nL 

with h the time-reversal breaking parameter. Assuming h ^ 0, we change variables to u = g/h and t = hhiL. Then 

Tt =1 + U ' (57) 

and direct integration gives 

u (t) = tan (t + arctanuo) ■ (58) 
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Thus the physics of the model repeats itself cyclically as the logarithm of the system size is changed. 
On the other hand, the functional conjugacy formalism gives 

*t (u) — exp ^ ™ 2 ^ = exp (arctanw) , "f^ 1 (u) = tan (In it) (59) 

u (t) = (A** (uo)) = tan (tin A + arctanw ) (60) 

Comparing the last expression to (|58|) we see that the Schroder eigenvalue for a unit step in t is A = e. 
In general, the step-scaling function corresponding to t step size A is 



a (u) = u (t)\ t=A = tan (A + arctanw) = — . (61) 

1 — u tan A 



„ . „ n , „ 1 + u da (u) 2 

(3(u) = l + u\ a(u) = - , = ~ -o, (62) 



u + tan A 

:tanw) = 

Choosing A = ir/4 for convenience, a quick check on (fTE)) gives 

1 +u 

1 - u ' du (1 - it) 2 
So j3 (a (u)) — da }^ P (u) indeed holds. In fact, the functional equation obeyed by ^, namely, 

a . , T / u + tanA A 

e A * (u = * r , 63) 

v 7 V 1 - u tan A / ' v ; 

actually belongs to the first class of examples discussed in §2 of the original paper by Schroder [23[ . (As an aside, 
we also note that ([ST|) is a complexified form of the Beverton-Holt-Skellam model from population dynamics flOl-) 

The trajectory (|6T)j) describes a limit cycle because it is periodic in t, as observed by LeClair et al. [17[ (also see [HI]). 
This is clearly a consequence of 'J -1 being periodic, as we remarked earlier in a general context, in §11. However, the 
discontinuity in the trajectory is somewhat peculiar to the model. This discontinuity is also displayed by the step 
scaling function, as shown in the Figure. 




-4 



' -5 " / 

Step-scaling function a (it) = tan (^ + arctannj, in blue, compared to the identity map. Note the discontinuity at 

u = 1 and the lack of a real fixed point. 

Since a does not intercept the identity map there is no real fixed point in this case, although ^ can be constructed 
by series solution about the purely imagin ary fixed points at u = ±i. 

We defer to the detailed discussions in [17| for a more complete picture of the physics described by this example. 
We return to a model-independent viewpoint, to explore other ways that limit cycles might be encountered. 
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V. MORE ON THE /3 FUNCTIONAL EQUATION 



For emphasis, we state again the functional equation relating the renormalization "velocity" j3 to the step-scaling 
function a. 

P(v(u)) = ^f3(u) . (18) 

Generally speaking, if j3 is known, this (nonlinear) functional equation determines a (u), up to a constant of integration, 
while if a (u) is known, the equation determines j3, again up to a (normalization) constant. However, the equation 
may hold in store some surprises. 

A fixed point in the step-scaling scheme must obey a (u*) = w*. Thus at a fixed point, it is trivially true 
that /3(cr(u*)) = /3(u*). But what has this to do with the standard fixed point condition that /?(«*) = 0? If 
da /du\ u ^ 1, then it would seem to follow from the functional equation that j3 (u*) = when a (u*) = u*, while if 
da /du\ u — 1, the functional equation itself would not lead to any conclusion about the value of /? (u*). But there 
is a subtle assumption here: All this is true if there is only one branch of the analytic function giving rise to f3. 
Related to this, a zero in da/du may induce, coincide with, or even supplant a zero of f3. This is at odds with the 
usual renormalization group point of view, but in general this too can occur, and when it does, it may foreshadow a 
much richer renormalization structure. 

If we have a zero in da/du, say at uq, and j3 (uq) is finite, then the functional equation implies (3 (a (uq)) = 0. 
Now, in the usual renormalization group situation, which we may describe as a purely first-order framework, u (t) is 
completely determined by u (to) and a single function (3 (u) = du/dt. So, if (3 (a (ito)) = and a (u) is just the initial 
u after having been i-evolved for some discrete step in t, this means that a (uq) = uq, and therefore (3 (uq) = must 
also hold. The fact that da / du\ UQ vanishes does not make much difference in this purely first-order point of view [3l| . 

However, in the quasi- Hamiltonian framework defined and discussed in Q, it may be that a (uq) is not just a zero 
of /3, but rather it is a branch point of an analytic function whose various branches constitute a family of (3s. In this 
situation, to completely determine the ^-evolution of u it is necessary to specify a transition function. That is to say, 
it is necessary to give a prescription describing how the trajectory switches from one branch of (3 to another when the 
branch point a (uq) is encountered. 

In this approach, a (uq) is a "turning point" in the evolution, and not necessarily a fixed point. Hence it is not 
necessary for a (uq) to be the same as uq. Nor is it necessary for (3 (uq) to vanish. 

The quasi-Hamiltonian approach brings to mind a Hamiltonian system whose underlying dynamics is actually 
a second-order differential equation, but which has been reduced to a first-order system between turning points 
through the use of energy conservation (see the first Appendix). Thus one would expect, at the very least, that the 
corresponding j3 function could flip sign at the turning point a (uq). 

In fact, the situation can be more complicated. There may be more than the two branch choices ± \/3\ at the 
turning point. In general there may be an infinite number of branches from which to choose j3\ 

A local, differential way to think about the various alternatives is to consider the RG acceleration, jerk, etc., along 
the trajectory as computed through use of the chain rule: 

^iu)_ Hu) d^-y_ { u)\ _ 



dt n v ' du V dt n 

So, for instance, if 



p{u) ~ (a(u )-u) p , (65) 



then 



— — — ~ (a(u )-uy w . (66) 
ar u->o-(uo) 

Thus, even though j3 (a (uq)) = 0, the RG acceleration d/3/dt can be nonzero at u = a (uq) if < p < 1/2, indicating 
that a (uo) is a turning point and not a fixed point for the continuous flow under these circumstances. 
All these things are best understood through consideration of some toy models. 
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VI. TOY MODELS AND NOVEL BEHAVIOR 



Consider a toy example that is unphysical (so far as we are aware) but can be solved in closed form [23], 

a(u) =2u(l-u) . (67) 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 

U 



Toy model a (u) = 2u (1 — u) (solid red) and a 1 (u) (dashed red) versus u. 

There are fixed points of a (u) at and = 1/2, as indicated by red dots in the Figure. The functional equation 
for the (3 function in this case is 

P {2u (1 - u)) = 2 (1 - 2u) P (u) , (68) 

with the normalization determined to be In 2 from the eigenvalue of Schroder's equation at u = 0, namely, 2^P (u) = 
'I' (2u (1 — u)). Series solution of the functional equation for /3 about u = immediately gives 

9 9 4 16 3? 1 
p(u) /ln2 = u-u 2 -=-u 3 --u 4 --u 5 ~— u 6 -?-u 7 + ■■■ = -(2u-l)ln(l-2u) , (69) 
3 3 5 15 21 2 

where the latter closed form is not difficult to guess from the first twenty or so terms in the explicit series, and is 
easily checked to solve the functional equation exactly. This /3 has zeroes precisely at the fixed points of a, again as 
indicated by red dots in the following Figure. 



beta(u) 




The toy P function, p (u) = \ (In 2) (2u - 1) In (1 - 2u). 
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For this special case, da/du = 2(1 — 2u) vanishes at it* = 1/2, the nontrivial fixed point of a, and in addition 
/3(l/2) = 0, but these coincidences are not true in general, as we shall see. In fact, in spite of the divergence in 
dj3{u) /du\ u _ u , both the RG acceleration and the RG jerk vanish at the nontrivial fixed point, (u*) = = 

4rzl3 (u*), as do all higher t derivatives of /3, since liniu^]^ {2u — 1) In™ (1 — 2u) — for any n. Thus in this case the 
RG flow into the nontrivial fixed point takes place with a very "soft landing." 
Integration of du/dt = (3 (u) gives 



u(t) = ±(l-(l-2u(0)f) 



This is precisely of the FC form sanctioned by Schroder, namely, u (t) 
and "f^ 1 (x) = | (1 — e x ). A representative trajectory is shown here. 



*- 1 (2**(m (0))), with 



(70) 
In (1 - 2x) 



0.6 - 

u(t) 




0.4- 




1 1 ! 1 1 


1 1 1 1 1 



-4 



-2 



A trajectory for the a (u) = 2u (1 — it) model, with u (0) = 1/4. 
More generally, consider toy models for step-scaling functions based on the logistic map 0, [ll|. For parameter s, 



with < s < 4, let 



er (u, s) = s u (1 — u) 



(71) 



This a has fixed points at and u* = 1 — - . Except for the special case s — 2, does not coincide with the location 
of the maximum of the map where da/du = 0, as given by a (1/2) = s/4. 
The functional equation for the (3 function in this case is 



/3 (su (1 - u) , s) = s (1 - 2u) {} (u, s) 



(72) 



where the normalization of the solution is given by ln(s). That is to say, (3 — (Ins) vp/ (d^f/du) where ^ is the 
corresponding solution of Schroder's equation with eigenvalue s. 



sf(u,s)=f (su(l-n), s) 



(73) 



While these functional equations do not admit closed-form solutions, for generic s, they can be solved numerically 
with sufficient accuracy for our purposes here through a combination of series and functional methods. 

One striking result of this numerical analysis is that (3 (u*) ^ 0, for 2 < s < 4, so the fixed point of the step-scaling 
function is actually not a true fixed point under continuous changes in scale! 

Explicit series solution of the functional equation for j3 about u = gives 



(3 (u, s) I (In s) = u 

4 + 5s + 7s 2 

> 7;' 

( s 2 _ !) ( a 2 + s + 1) ( s 2 + !) 

8 + 10s + 21s 2 



1 



4 + 5s 



1 

2- 



1 



(s 2 - 1) (s 2 + s + 1) 



18s + 31s 2 + 42s 3 + 35s 4 + 21s 5 



(s 2 - 1) (s 2 + 1) (s 4 + s 3 + s 2 + s + 1) (s 2 + s + 1) 
25s 3 + 39s 4 + 21s 5 + 33s 6 



(s 2 - 1) (s 2 + 1) (s 4 + s 3 + s 2 + s + 1) (s 2 + s + 1) (s 2 - s + 1) 



Oi 



(74) 
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More generally we write 



P (u, s) I (In s) = u ^1 + ^ Cn (s) u' 1 ^ , Ci = 



C-2 = 



1-S 2 ' 



(75) 



with higher coefficients in the series given by the recursion relation (here [""J is the integer- valued floor function) 



C n +2 (s) 



I — s n+2 



^ " +1 / 1 \ \ 

2c n+1 (,s) + L/-1- K-^^W forn>l. 



(76) 



Numerical study § of c„ (s), for various values of s and n < 200, provides compelling evidence that ([75]) converges 
for \u\ < R(s) where 



R(s) 



lim sup (|c„ (s)| 1/n 

n— > oo \ 



if < s < | , 



if I < a < 2 , 



(77) 



if 2 < s < 4 



A closed form is not known for (|75jl. except for s = 0, ±2, and 4. Nonetheless, as already mentioned, the model is 
amenable to numerical analysis for generic s. 

In particular, the functional equation can be exploited to continue the series and exhibit the various branches, j3 n , 
of the multi- valued (3 function that is encountered for 2 < s < 4 0. For example, 



Po (it, s) — \/ s 2 — 4su [3 ^— — V s 2 — Asuj , s^j , 
Pi (u, s) = -V ' s 2 - 4su /3 ^ (s + ^s 2 - 4su^) , s^j 



(78a) 
(78b) 



where /3 is the explicit series (|75[) , /3q is the continuation of this series through use of the functional equation to give 
the principal branch on the interval, < u < s/4, and /?i is the first alternative branch which is real- valued on the 
sub-interval js 2 (l — is) < u < 4. Etc. For s < 2 only one branch is needed, namely, /3o, but additional branches, 
such as pi, are required to develop completely the trajectories for s > 2. 
of real- valued branch functions is given by iterating the definition of Pi . 



In the latter situation, an infinite sequence 
Thus, 



P n+ i (u, s) = -V 's 2 - 4su p n f i (s + V 's 2 - 4suj , s^j 



(79) 



These are all the additional branches of P that are needed for 2 < s < 3 and for s = 4, but for 3 < s < 4 there are 
other branch function sequences that must be taken into account to describe fully the continuous trajectory u (t) . 

By construction, stemming from the fact that a is quadratic, for s > 2 the branches of P given by (I79| have 
square-root zeroes at points obtained by iterating the action of the step-scaling function, starting from uq = 1/2. 
That is, zeroes are given by the sequence {a (uq) , a (a (ito)) ,•••}• Thus the p n branches arrive at those zeroes with 
infinite slope, exactly as described by (|65[) . and its u derivative, with p = 1/2. On the other hand, the product 
P dp/du = dp/dt is finite and nonvanishing in the limit as u goes to one of these zeroes. Consequently, the RG 
acceleration does not vanish for points in the sequence of P zeroes: They are turning points, not fixed points. 

Further general discussion of this class of toy models would take us too far afield. Suffice it to consider here three 
other explicit examples, two based on numerical analysis (s — 11/4 and s — 10/3), and one based on elementary closed- 
form expressions (s = 4). For convenience, we first take s = 11/4. This gives a step-scaling model a (it) = -M- u (1 — u), 
with fixed points at and w* = 7/11. 

Constructing graphs like those for the previous toy model, we find that there are several interesting differences 
between the Figures for the two models. The slope of the step-scaling function is now negative at the nontrivial fixed 
point, as opposed to the vanishing slope of the previous toy model. This negative slope is a crucial ingredient that 
gives rise to the multi- valued-ness of the P function, as shown in the Figures. 
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0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

U 



Toy model a (u) = ^ u (1 — u) (solid red) and a 1 (it) (dashed red) versus u. 




-0.2 -- 



Six branches of the multi- valued (3 function for the model with a (u) — 4- u (1 — it). 

To avoid confusion, we emphasize that this last "phase-space" Figure shows the real- valued branches of /3(u), i.e. 
the Figure does not show a trajectory in a two-dimensional coupling space. The trajectory u (t) is one-dimensional, 
just as it was in the previous toy model, as illustrated in the following graph. Although u* = 7/11 is a fixed point of a, 
it is not a zero of /3 on any of its real- valued branches. Moreover, at the nontrivial zeroes of /? the RG acceleration does 
not vanish because the relevant branch of (3 has a compensating infinite slope, as discussed more generally following 
([79]) . and as is evident in the plot of the toy trajectory, so these zeroes are indeed turning points. Nevertheless, the 
net effect of the multiple branches of /3 on the trajectory u (t) is to produce an oscillatory convergence to it* = 7/11 
as t —¥ oo. This is distinctly different from the monotonic approach to u* = 1/2 of the previous toy model whose f3 
function had only one real branch. 
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-4 -2 2 4 

t 

A trajectory with u(0) = 1/4, for the s = 11/4 multi- valued ft function. 

Next, consider s = 10/3, another case which yields to numerical analysis but cannot be described with closed- form 
results. The step-scaling function in this case is a (u) — u (1 — u), with fixed points at and u* = 7/10. The 
graph of cr is similar to the previous s = 11/4 example, and is left to the reader to plot. This toy example is especially 
interesting in that it provides an explicit RG realization of a limit cycle, a possibility conjectured by Wilson in his 
early, classic study 26], but not thought to be possible for a model with only one coupling until relatively recently 
[27j . Here, the exact cycle is only realized asymptotically in the limit of very large t. 

The basic mechanism whereby this is achieved for the s = 10/3 case is the same as in the previous toy model: There 
are an infinite number of branches of the underlying analytic ft function, and the evolving RG trajectory switches 
from one branch to another when turning points are encountered. The branch structure is more complicated for 
s = 10/3 than for s = 11/4, however, with the branch points accumulating around ui ow = || — ^j\/l3 = 0.4697 and 
Mhigh = ^j\/l3 + |§ = 0.8303 in such a way that, as t — > 00, the trajectory approaches a rectangular sequence of 
steps between these two values, hence giving rise to a two-cycle asymptotically in t. The initial stages of this large t 
behavior are evident in the following sample trajectory. 



U(t) 0.8 - 






r 


0.6- 








0.4 - 








< 1 1 1 1 


1 1 1 1 1 1 1 



-4 -2 2 4 6 

t 



A trajectory with u (0) = 1/4, for the s — 10/3 multi- valued ft function with a 2-cycle limit. 

While the methods used in the previous example are also applicable to this case (see (|78a[) . (|78b[) . (|79|) . and 
generalizations), we will not construct here the actual branches of ft for s = 10/3, nor will we enumerate the sequencing 
of the various branches as the trajectory evolves. This information can be found, in a different context, in Q. Suffice 
it to say that although u* = 7/10 is a fixed point of cr, once again it is not a zero of ft on any of its real-valued 
branches. 
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For a final, chaotic example, which nevertheless admits an exact, closed- form solution (23|, consider 

a (u) = Au (1 — u) , 

for < u < 1. In this case the real- valued branches of the analytic /3 function are given by 

1 + n 



/3„ («) = (In 4) y/u(l-u) ((-iy 



n + arcsm \/u 



for all integer n > 0. The functional equation is indeed obeyed, in the following sense: 

f) n (cr (u)) = p n (u) da (u) /du for < u < 1/2 , 
/3 n+1 (cr (u)) = j3 n (it) dcr (u) /du for 1/2 < u < 1 . 



(80) 



(81) 



(82a) 
(82b) 



Note that cr (u) approaches the branch point at cr = 1 as u goes to 1/2. Beyond this, for 1/2 < u < 1, the j3 (a (u)) 
term in (p~8|) switches from the nth branch to the (n + l)st. Correspondingly, at the nontrivial fixed point = 3/4 
of the step-scaling function, we see that j3 (m*) does not vanish for any of the branches of f3. A typical RG trajectory 
in this chaotic case is shown here. 




A trajectory with u (0) = 1/4, for the s — 4 multi- valued /? function leading to chaotic evolution. 

There is one highly unusual feature for the trajectory shown, as well as for all other trajectories for this example: 
The coupling actually goes to zero with increasing frequency but always "bounces back" to positive values. For 
a thorough discussion of this model using functional methods, in the context of classical mechanics for a chaotic 
dynamical system, see @. 



VII. DISCUSSION 



In this paper we have elucidated, explored, and applied the underlying functional conjugacy structure of the 
renormalization group implicit in Gell-Mann and Low's finite renormalization group equation [14j |. structure which 
is normally overshadowed by local differential (Lie algebraic) features. We introduced methods to extract the links 
among fixed points of exemplary RG trajectories in section II. We applied some of these functional methods to obtain 
continuous flows from step-scaling functions utilized in lattice gauge theory, in section III. We explained in section 
IV how novel features could arise in general. We illustrated such features — including multi-valued /? functions, limit 
cycles, and chaotic trajectories — using toy models based on the logistic equation in section V. 

The functional conjugacy core of the RG, ([6]), amounts to a solution by the method of characteristics, (|T2|) ct 
seq., so that scale changes are equivalent to variations of the initial coupling data. This controlling characteristic 
structure implies nonlocal associations in renormalization flows. In the large, the resulting functionally conjugate 
RG trajectories possess an elegant global mathematical consistency which sheds light on the flow of couplings in the 
presence of both UV and IR fixed points, and which also shows that trajectories may admit turning points, i.e. zeros 
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of (3 functions where not all t derivatives of (3 vanish as a consequence of intricate branch structure. This latter 
possibility may lead to limit cycle or even chaotic behavior, as we have illustrated in some detail. 

In this paper, models with only a single RG coupling flow have been considered . Extensions of the functional 
methods to models with more than one flowing coupling have yet to be carried out in complete detail. This is an 
area that warrants further exploration. 

While we have not yet fully explored realistic quantum field theoretical models evincing all of the possibilities 
that we have discussed, we are confident that they do exist, at least for many of the features we have described. 
Presumably, the more exotic, oscillatory RG flows would only emerge in systems where the c-theorem [a, l28j (or its 
equivalent in higher dimensions [2l[) does not hold [ill H3- We would be most pleased to find in actual physical 
systems the full range of behavior illustrated by the toy models. 
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VIII. APPENDIX: NEWTONIAN TRAJECTORY AS A TRANSPORT OF DATA 



Here we discuss analogies between RG flow and one-dimensional motion of a classical particle, thereby making 
contact with our previous work 0-0 • 

For fixed energy, E, the trajectory x(t, Xq) of a particle moving in a one-dimensional potential, between turning 
points, depends only on time t and initial position xq — x\ t=0 , since the initial momentum is fixed by E up to a choice 
of branch for y/E — V (x) (usually just an overall ± sign if V itself has only one branch). Thus 

= t, (83) 



v (x) 

where v (x) is the velocity profile along the trajectory. In this context it is clear that v = is not necessarily a fixed 
point. Rather more often it is a turning point of the motion. 

Now, here's a simple technique that appears in Gell-Mann and Low (see Eq (B20), [Hj]) : Differentiate ([83]) with 
respect to the initial position, regarding t and x$ as independent variables. This gives 

1 dx(t,x ) 1 _ Q 



v (x (t, x )) dx v(x ) 
On the other hand, v (x (t, xq)) — dx (t, xq) /dt, so this last result is just 



dx(t,x ) dx(t,x ) 

-v{x ) — . (85) 



dt dxo 

That is to say, for want of a better name, this is a one-dimensional "Gell-Mann-Low transport equation" for x (t, xq 

JLf — JL 

dt J dx 



It is not quite "advection" |3C| since that would have the form of a conservation law, namely, -Srf — -§- (vf). The 



two types of transport are exactly the same only for constant v. 

In fact, (1851) is simpler than advection, with solutions of the FC form 



x(t,x ) = y- 1 {e t y(x )) , (86) 

for an appropriately defined Schroder function ^ (although this terminology is not used in [3) with inverse function 
\f r_1 . By "appropriately defined" we mean, 

«M = *to)/*M = pM^y^y (^) 
hence \f r is essentially just the exponentiated time to reach x from some reference point, 

y(x) = *(x Icl ) CX p(J* A) . (88) 
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Another way to express the result (|86l) is as a formal Taylor series in t, rewritten in terms of v (xq) and its derivatives 
through the use of (|85|) . Thus 



°° I gn 00 1 / Q \ 



rw <9t 

n=0 



71 

r=0 ra=0 



1 / 5\ 71 ^ 

t v (x ) + E ~i *" r ( x °) JT~ v ^ ■ (89) 



n— 2 x 

(See Eq (B21) in [3.) If it is not already evident why this solution is formally the same as then repeat some 
steps from the text and write the series in (|89[) as exponentiated operators. 

x(t,x ) = x(t,x ) = e v{xo) a; = e* 81 »*(-o) a; , (90) 

where in the last step we have used (|B7|) . But now, a; = (* (zo)) = ( ex P (In* (#o)))> so the last expression 
in (|9"0|) reduces to a translation of the variable In (* (&o))- 

e * am^o, $-1 ( ex p(ln*(x ))) = (exp (t + ln*(a:o))) = * _1 (e**(a; )) . (91) 

Thus we recover (|86|) as an element of an abelian Lie group, given as usual by exponentiating an element of the 
underlying algebra — a symplectomorphism in the present context of classical mechanics. 

IX. APPENDIX: A VARIABLE CHANGE AT A LATTICE FIXED POINT 

For a supposed nontrivial fixed point, u* = 1/gl, determined by 

er (u*) = u* , (92) 
with f3 («*) = 0, we change variables in the lattice data formulas (|42j) and (j43]) of the text, by writing 

- s =(l + r)gl, (93) 
so that r > for bare couplings <?q above the fixed point value g^. Under this change of parameterization, we have 
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(£) (1 + rY g^i = «.+$> (I) r j , (94) 
j>i J 

= £^> C4 = (i + l)a i+1 g, | = -^. (96) 



ax = 4- -l-^0'-l)q, WJo 2 j I . 02 = 4- (l-^E^- 1 )^- 2 )^'^)^) , (97) 

»kmax(n) = 4" ( (-1)" ~ ^ E 0' " 1) U ~ 2 ) ' ' ' U ~ *) 9 (*) So* ) , 0* = . ( 9 «) 



and therefore p- Sjl n +i ( — 1) J rJ = (i+r)g^ i as it should. Similarly 



<r («)=«♦+$> ^^ = E a H i ) ? - i - (99) 



j>l j>0 
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Discarding an overall 4^, the functional equation for the (3 function (|44l) is then 

]T U + 1) <h+i W rJ J Z 3 ( "* + E Qfe ( L ) r * I = f E C7 + 1) W r 1 ' ) /3 U* + E n * w r " I • ( 10 °) 
V>>° / \ fc>i / \i>o / \ k>i J 

Now it is straightforward albeit tedious to construct, up to an overall normalization, a series solution in r of this 
functional equation to obtain an expression for j3 near the fixed point. We shall assume u* is a first-order zero of fi. 

min(n, N) 

p= E b * r * ■ ( 101 ) 

i=l 

It would be nothing more than wishful thinking to carry the series approximation for j3 beyond that for either u or 
a(u). 

For numerical work described in the text, we build the series solution to fourth order in r. Define a.j (£) = ttj, 
dj (L) = and 

/3 (u* + w) = biw + b 2 w 2 + b 3 w 3 + b 4 w i . (102) 

Then (fTOOj) gives 

(oi + 2o 2 r + 3a 3 r 2 + 4a 4 r 3 + 5a 5 r 4 ) /? (u„ + Stir + Sl 2 r 2 + Sl 3 r 3 + St 4 r 4 ) 

= (2li + 2St 2 r + 3St 3 r 2 + 4Sl 4 r 3 + 53t 5 r 4 ) (3 (u, + a x r + a 2 r 2 + a 3 r 3 + a 4 r 4 ) . (103) 

Now, we obviously cannot determine the overall normalization of j3 from the functional equation (|100j) . This fact is 
why we find the same coefficient of r 1 on the LHS and RHS of this last equation, for any b\. But higher powers of r 
give nontrivial information. For example, r 2 gives 



b 2 _ 2lia 2 - Stall 
b~i ~ Slioi (a 4 -2li) 

r J gives 



(104) 



bs Sl 2 Q 2 -Sl 2 a 2 + St 1 a 1 (gi_a 3 - Sl 3 Oi 
bi St 2 a 2 (a 2 - St 2 ) 



2 o,2„2 /„2 ol2^ > ( 105 ) 



and finally, r 4 gives the unwieldy expression 

2Sl?St 2 aiai + 2Sl 2 Sl 3 a?a 2 - 2St?St 2 a 2 a 3 - 2StiSt^a?a 2 - 8St?a 2 o 2 a 3 + 8Sl 2 Sl 2 Sl 3 a? 
-7Stfaio 2 o 3 + 7StiSt 2 St 3 af + 3St 4 a 2 a 4 + 3St?a?a 4 - 3St?St 4 a? - 3Sl 2 St 4 a 4 
+5St?aiai - 5StiSt|af + 4Stfo| - 42t|af + St?St 3 a 2 a 2 - St 2 St 2 a?a 3 - St 2 St^o 2 a 2 + Sl 2 St 2 a 2 a| 

Sl 3 a 3 (a 2 -Sl 2 )(a 2 + a 1 St 1 +Sl 2 ) 



(106) 



And so it goes. We obtain bj/bi as functions of the as and Sis. But again, it is overly ambitious to carry the series 
approximation for /3 beyond that for either u or a (u). 

There remains only one coefficient to be determined, namely, bi. This is given by 

bl=mA = m f^P) . (107) 



oi (£) 

Although repetitious, a complete derivation of this result goes as follows. Setting the scale of t so that u(t = 1) = a (u), 
Schroder's equation is, once again, 

#(<r(u)) = A*(u) , (108) 

and from u (t) = f" 1 (A** (u)) with ^ = P (u (t)), we obtain /3 (u (t)) = (In A) * (u (t)) (u (t)). Now, at a fixed 
point u* = a (u*), with VP (u») = 0, a series solution of (|108j) in powers of (u — u*) requires 



da (u) 
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du 



da/dr 



du/dr 



r=0 



ai ( £ ) (109) 
oi (£) ' (iUJj 
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and also that j5 (u) = (In A) (u — u*) + ((u — it,) 2 j . Hence (|107[) is obtained. Alternatively, in terms of the original 
series involving the bare lattice coupling, we have 

x _ 1 + Ej>2 (j - 1) c j ( L ) go* _ u * + Ej>i i c i ( L ) floT 1 , > 

~ 1 + E 3 >2 (i " 1) 9 W 3o J * " u. + E,>! J 9 (I) git 1 ' 
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